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A Remark on Formation Control with Triangulated Laman Graphs: 
Genericity of Equivariant Morse Functions 
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Abstract —This paper, as a continuing work of [1], focus on 
establishing the fact that if we equip a reciprocal multi-agent 
(RMA) system with a triangulated Laman graph (TLG), then 
the associated potential function is generically an equivariant 
Morse function, i.e., there are only finitely many critical orbits 
each of which is nondegenerate. Though this assumption on the 
potential function of being an equivariant Morse function has 
been used, and in fact indispensable, in several occasions. But 
it is actually still an open question whether it is true for a given 
RMA system. Thus, in this paper we will provide a confirmative 
answer to the question for the class of RMA systems with TLGs. 
The main result, as well as the analysis of this paper, has many 
implications for other difficult problems. 

I. Introduction 

In this paper, we will consider a special class of reciprocal 
multi-agent (RMA) systems equipped with a particular type 
of Laman graphs, as we call the triangulated Laman graphs. 
It is known that each RMA system is a gradient system with 
respect to an equivariant potential function, and in this paper 
we will develop among other things a basis property of this 
special class of RMA systems. We show that for each of such 
systems, the associated potential function is generically an 
equivariant Morse function, i.e, there are only finitely many 
critical orbits each of which is nondegenerate. 

We start by reviewing the mathematics model of a RMA 
system. Let < G = (V 1 E) be an undirected graph of N vertices, 
with V = {l,--- ,A'} being the set of vertices and E the set 
of edges. We say two vertices i and j are adjacent if (7, j) 
is an edge of G, and for convenience we let V? be the set of 
vertices which are adjacent to vertex i. Let {/,-,j (7, j) £ £} be 
a family of continuous differentiable functions from R + , the 
set of positive real numbers, to M. The equations of motion 
for a set of N agents jci, - - - ,xn are described by 

Xi = E fiMij) ■ ( Xj-Xi ), Vi = 1, ■ ■ ■ ,N (1) 

JeVj 

where djj denotes the Euclidean distance between x, and Xj. 
As we see from the equation above, the interaction between 
Xj and x; is modeled by and it depends only on the 
distance. Along this paper, we assume that fjj = /),- for all 
(7,7) £ E, i.e, interactions among agents are reciprocal. 

It is known that each RMA system is a gradient system, 
with the potential function given by 

rdi j 

,X N ) := £ / xfjj(x)dx (2) 
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It is a summation of pairwise potentials between adjacent 
agents. Notice that the potential function <b depends only on 
relative distances among agents, so the value of <f> is invariant 
if we translate and/or rotate the whole configuration. In 
mathematics, it says that <t> is an equivariant function with 
respect to the group action of rigid motion. A precise 
definition will be given later in section III. 

An open question left so far is to ask whether it is 
true, or at least generically true, that the potential function 
<t> is an equivariant Morse function, i.e, there are only 
finitely many critical orbits (orbits of equilibria) each of 
which is nondegenerate? This question is important not 
only because the mathematical statement itself, but also 
because it is an indispensable condition assumed in many 
developed theorems and their applications. For example, 
analyzing local exponential stabilities, counting number of 
critical orbits [2], characterizing their regions of attractions, 
investigating system behaviors under perturbations [3]—[5] 
and etc., all these questions rely on the assumption that <t> 
is an equivariant Morse function. 

In this paper, we will prove that if the underlying graph 
G is a triangulated Laman graph (TLG), as we will define 
in the section III, then generically <f> will be an equivariant 
Morse function. We have showed in [1] that a RMA system, 
equipped with a TLG, has several distinct properties. For 
example, for each equilibrium configuration of system 0. 
there is a geometric decomposition of the configuration 
into union of line sub-configurations each of which is also 
an equilibrium. Furthermore, we have developed a formula 
which relates the Morse-Bott index of a critical orbit, as 
an algebraic term, to this geometric decomposition. In par¬ 
ticular, the formula says that the Morse-Bott index of a 
critical orbit can be computed as the sum of the Morse- 
Bott indices of critical orbits of these decomposed line sub- 
configurations. This formula has a potential impact on the 
design and control of RMA systems as it enables us to locate 
or place critical orbits with various Morse-Bott indices over 
the configuration space, we here refer readers to [1] where we 
apply this formula to count and locate all the stable critical 
orbits within a special class of interaction laws. 

In this paper, we will follow the results achieved in paper 
[1] to develop several other relevant properties associated 
with this special class of RMA systems. After this introduc¬ 
tion, we proceed as follows. In section II, we will introduce 
a class of interaction laws which will be considered in 
this paper, and also we will introduce the Whitney C 1 - 
topology on it. In section III, we will review some key 


definitions about RMA systems, and then we will state the 
main theorem of this paper. This theorem claims that the 
potential function <t> is generically an equivariant Morse 
function. The rest of the paper is then devoted to the proof of 
this fact. In section IV, we will briefly review the canonical 
partition and the Morse-Bott index formula, as two key 
results developed in [1]. In section V through VIII, we will 
mainly focus on critical line configurations. In particular, we 
will show in section VIII that generically, there are only 
finitely many critical orbits of line configurations each of 
which is nondegenerate. By combining this result with the 
Morse-Bott index formula, we will then be able to establish 
the genericity result. 

II. Class of monotone interaction laws 

In this section, we will introduce a class of interaction 
laws, together with the C 1 -Whitney topology on it. 

Let G be an undirected graph of N vertices. Let P be the 
configuration space defined by 

P:= {(jfi,— ,x N ) €R 2xN \xi^Xj,W(i,j) £E} (3) 

Configurations with collisions of adjacent agents are ex¬ 
cluded. We will now define a class of interaction laws by 
which the solution of system Q. with any initial condition 
in P, exists for all time, and converges to the set of equilibria 
of system (ID- 

Let R be the set of positive real numbers, and let 
C (R +; R) be the set of continuous differentiable functions 
from R + to R. For each function / in C^R+jR), we let 

f(d):=df(d) (4) 

We introduce f because if / is an interaction law between 
a pair of adjacent agents, then / represents the actual 
magnitude of attraction/repulsion between them. We will use 
/ and / in various occasions, and both of them are useful 
in this paper. 

Let & be a subset of C 1 (R + . R) defined as follows. A 
function / is in & if and only if 

Cl. f(d) > 0 for all d > 0, and / has a (unique) zero. 

C2. lim [ f(x)dx = <=°. 
d-^oJd 

We here note that the two functions / and f share the same 
zero, i.e, f(d) = 0 if and only if f(d) = 0. We impose these 
two conditions because the first condition implies that the 
interaction is a repulsion at a short distance, and an attraction 
at a long distance. The second condition prevents collisions 
of adjacent agents along the evolution, so then the solution of 
system 0- with any initial condition in P, exists for all time. 
Moreover, we have showed in [6] that if each interaction law 
fij satisfies conditions Cl and C2, then all critical orbits of 
system (JTJ are contained in a compact subset of P. In this 
paper, we will assume that each interaction law fij is in 
thus we assume the convergence of system 0- 

We will now equip this function space & with the Whitney 
C 1 -topology. First we will describe the Whitney C 1 -topology 


on C*(R+,R), and we describe it by defining a basis of 
open sets in C*(R+,R). Let C°(R+) be the set of continuous 
functions from R + to R + . Then for each / in C^R+jR), 
and each 8 in C°(R+), we define an open ball Bg(f) of 
C'(R+,R) by collecting any g in C ! (R+,R) such that 

\ g (d)-m\+\ g '(d)-f\d)\<8(d) ( 5 ) 

for all d > 0. By varying / over C*(R+,R) and 8 over 
C°(R+), we then get a basis of open sets in C^R+jR). 

The subset C C 1 (R + ,R) is then equipped with the 
subspace topology. Notice that & is an open subset of 
C 1 (R+,R) with respect to the Whitney C 1 -topology. So a 
subset V of & is open in jF if and only if it is open in 
C'(R+,R). 

Let |£j be the cardinality of the edge set E, and in the 
case where G is a Laman graph, we have \E\ =2N — 3. 
The ensemble of interaction laws of system 0 should be 
considered as an element in In this paper, we equip 

with the product topology, i.e, a set U of is open 
if and only if it is a union of FT { j\ eE Vij with each Vij an 
open set in &. Finally, we say a property is generic if it 
holds on an open dense subset of 

III. Definitions and the main theorem 

In this section, we will review some definitions, as well 
as notations introduce in [1], and we state the main theorem 
of this paper. 

1. Triangulated Laman graph. Laman graphs are knowns 
as the minimally rigid graphs in R 2 [7], and each Laman 
graph can be constructed via a Henneberg construction. A 
triangulated Laman graph (TLG), as we will define now, can 
be constructed via a special Henneberg construction. Start 
with an edge, we then join a new vertex, at each step, to two 
adjacent existing vertices via two new edges. An example of 
a TLG is illustrated in figure [T] 



Fig. 1. An example of a TLG. Start with edge (1,2), we then subsequently 
join vertices 3, 4 and 5 to two existing adjacent vertices. 

2. Group action of rigid motion. Let SE( 2) be the special 
Euclidean group for R 2 , each element y in SE( 2) can be 
represented by a pair (0,v) with 0 in the special orthogonal 
group 519(2) and v a vector in R 2 . In this representation, 
the group multiplication of two elements yi = (0i,vi) and 
y 2 = (0 2 ,V2) is given by y 2 • yi = (020t,02Vt+V 2 ). 







We now define a 5’£’(2)-action on P by sending y = (0 .v) 
in SE (2) and p = (x\ , ■ • • ,xp) in P to 

y-P'= (9xi+v,--- ,0x N + v) (6) 

This group action is often referred as the group action of rigid 
motion because it preserves the shape of a configuration. In 
this paper, we let 

0 p -.= SE{2)- P (7) 

be the orbit of p with respect to the 5’£’(2)-action. 

3. Equivariant Morse function. An important observation of 
the potential function <t> is that <1> depends only on relative 
distances between agents. Consequently we have 

&(P)=*(Y-P) (8) 


Inductive step. Suppose G' — (V',E') is a subgraph of G 
consisting of vertices V' = {1, • ■ • ,n — 1}, and we have 
partitioned E' into disjoint subsets as 

E' =E[u---UE' ml (10) 

Now suppose vertex n joins to vertices i and j via edges 
(i.n) and and we describe the rule of updating the 

partition by taking into account (i. n) and (j,n). 

Without loss of generality, we assume that the edge (i.j) 
lies in E[, then there are two cases: 

Case I. If x u x) and x. n are aligned, then we update the 
partition by adding (i,n) and (j,n) into E[. 

Case II. If Xi, Xj and x n are not aligned, then we update the 
partition as follows 


for any p € P and any ys SE(2). In particular, if p is an 
equilibrium of system 0 - then so is p' in <& p . In any of such 
case, we call & p a critical orbit. 


Let f/ p be a critical orbit, and let It p be the Hessian matrix 
of 4> at p, i.e. 




d 2 <S>{p) 
dp 2 


(9) 


Since interactions among agents are reciprocal, the Hessian 
matrix H p is thus symmetric, and hence all eigenvalues of II p 
are real. The null space of H p at least contains the tangent 
space of Cj p at p. In other words, if we let no(H p ) be the 
number of zero eigenvalues of H p , then no{H p ) is at least 
three. We say a critical orbit f/ p is nondegenerate if no(H p ) 
is exactly three. (As the set of eigenvalues of H p i is invariant 
as p' varies over so this definition doesn’t depend on the 
choice of p.) 

A potential function <t> is said to be an equivariant Morse 
function if there are only finitely many critical orbits of 
system 0 each of which is nondegenerate. We will now 
state the main theorem of this paper. 


Theorem 1: Let G = (V,E) be a TLG, and we assume that 
the ensemble of interaction laws of system 0 is contained 
in The potential function 4> defined by equation |2| is 
then generically an equivariant Morse function. 


IV. The canonical partition and 

THE MORSE-BOTT INDEX FORMULA 

In this section, we will review some key results developed 
in [1], 


A. The canonical partition 

Let G = (V.E) be a TLG, and let p be a configuration 
in P. We will now introduce the canonical partition of E 
associated with p. Choose a Henneberg construction of G, 
and we label the vertices with respect to the order of the 
construction. The partition is then defined inductively by 
following the construction. 

Base case. Start with the subgraph G' = (V',E') of G 
consisting of vertices V' = {1,2}. Since there is only one 
edge (1,2) in E', the partition of E' is trivial. 


E[c • • ■ CE r m i U {(;,«)} U {(j,n)j (11) 

By following the Henneberg construction, we then derive the 
canonical partition of E associated with p. 

We proved in [1] that this partition does not depend on 
the choice of the Henneberg construction. An example of the 
canonical partition is illustrated in figure [2] 



Fig. 2. An example of the canonical partition. We see from the left 
figure that the graph G is a TLG as we label the vertices with respect 
to a Henneberg construction, and p is a planar configuration with xi,* 2>*3 
aligned, and X 3 ,* 4,*5 aligned. Then the canonical partition of E associated 
with p is given by the right figure. 


We will now list some relevant properties of the canonical 
partition. Let E = £j U • • • UE m be the canonical partition 
associated with p. Let G; = ( Vi,Ej ) be the subgraph of G by 
restricting G to £), and let p, be the sub-configuration of p 
associated with G;, then 

a) . Each G; is a TLG. 

b) . Each pi is a line configuration. 

c) . If there is another partition of E satisfying conditions a) 
and b), then it is a refinement of the canonical partition. 
In other words, the canonical partition produces minimal 
number of sub-graphs regarding to the first two conditions. 

d) . if p is an equilibrium, then each p,- is an equilibrium of 
the sub-system induced by G;. 

More details, including the proofs of these statements, can 
be found in [1], 





B. The Morse-Bott index formula 

In this part, we will follow the canonical partition to 
introduce a formula which can be used to compute the 
Morse-Bott index of a critical orbit. We start by introducing 
some necessary definitions. 

Let i & p be a critical orbit, and let H p be the Hessian matrix 
of <f> at p. Let n + (H p ), no(H p ), and ti-(H p ) be the numbers 
of positive, zero, and negative eigenvalues of H p respectively. 
In this paper, we refer to the triplet (, n+(M),no{M),n-(M)) 
as the inertia of H p . The Morse-Bott index and co-index of 
G p are defined to be n-(H p ) and n+{H p ) respectively. 

The Morse-Bott index formula. Let G be a TLG, and let 
G p be a critical orbit of system ([TJ. Let {pi}”T\ and {G,}'" = | 
be sub-configurations of p and subgraphs of G respectively, 
associated with the canonical partition. Let <t>, be the induced 
potential function associated with G,-, i.e, 

’S'lip'i) '■= L [ 1 fjk{x)dx ( 12 ) 

U,k)eEi Jl 

and let H Pj be the Hessian matrix of <t>,- at p,. Then we have 

( n-(H p )=ZT =1 n-(H Pi ) 

\ n + (H p )=Lr =] n+(H Pi ) 

This set of equations will be referred as the Morse-Bott index 
formula. This formula has some relevant implications as we 
state below 

1) . The critical orbit G p is nondegenerate if and only if 

each G Pi is nondegenerate. 

2) . Suppose G p is nondegenerate, then G p is (exponen¬ 

tially) stable, i.e, n + (H p ) =0 if and only if each G p , is 
(exponentially) stable. 

We refer readers to [1] for a complete proof of the formula, 
as well as the proofs of the implications above. 


V. The Hessian matrix at a line configuration 

As we have seen in the last section that the Morse- 
Bott index of G p can be computed as the sum of Morse- 
Bott indices of G P1 ,--- ,G Pm with pi,-- - ,p m the line sub¬ 
configurations of p associated with the canonical partition. 
So in this section, we will focus on the case where p itself is 
a critical line configuration, and we will compute the Hessian 
of 4> at p. 

In the rest of this paper, we assume that the a-axis and 
the b- axis are the two axes of R 2 , and we let a,- and bj be 
the two coordinates of 3c,-. Let 

( a:= (aw- ,a N ) 

\ b:= (, b \,-■• ,b N ) 

We then re-arrange entries of a configuration p so that 

p = (a,b) (15) 

The Hessian matrix H p will then be computed with respect 
to this order. 


Since the only information we need in this paper is the 
set of eigenvalues of H p , and since the set of eigenvalues of 
H p / is invariant as p' varies over the orbit G p , we may rotate 
and/or translate p if necessary so that the entire configuration 
p is on the a-axis. 

To compute H p , we will now introduce two N-by-N matri¬ 
ces F p and dF p . Both of these two matrices are symmetric, 
of zero-row/column-sum. So we define these two matrices 
by specifying their off-diagonal entries. Let F p l j be the i /-th 
entry of F p , and dF pjJ the ij -th entry of dF p , we then define 


and 



if (ij) e E 
otherwise 


dFp,ij — 


f'ij(dij ) if (i,j) S E 

0 otherwise 


The Hessian of <t> at p is then given by 



(16) 


(17) 


(18) 


We will now discuss about the null space of H p . The null 
space of H p at least contains T p G p , i.e, the tangent space of 
G p at p, as we will compute now. 

Let e be a vector in R. N of all ones, and we define two 
vectors in R 2xN by 


t a := (e, 0 ) 
t b := ( 0 ,e) 


(19) 


These two vectors represent the infinitesimal motions of 
translation along a-axis and ft-axis respectively. We now 
define another vector in Mr xN by 


r p := ( 0 ,a) 


( 20 ) 


This vector then represents the infinitesimal motion of clock¬ 
wise rotation of p around the origin. 

These three vectors t a , Tb and r p form a basis of the tangent 
space T p G p . By direct computation, we verify that all the 
three vectors t a , Tb and r p are in the null space of H p . On 
the other hand, if the critical orbit G p is nondegenerate, then 
the null space of H p coincides with T p G p . In particular, we 
note that if G p is nondegenerate, then the null space of dF p 
should only be spanned by e, and the null space of F p should 
only be spanned by e and a. 


VI. The reduced system 
In this section, we will introduce the notion of reduced 
system which will be a useful tool for analyzing and com¬ 
puting the inertia of the Hessian matrix at a critical line 
configuration. 

We start by introducing the notion of virtual interaction. 
Consider an auxiliary system of three agents 3ti, X 2 and X 3 . 
Let /12 and /13 be interaction laws between x\ and j? 2 , and 
between x\ and X 3 respectively. We assume that xi, X 2 and X 3 
are aligned on the a-axis, and we will now define the virtual 
interaction between X 2 and X 3 induced by /12 and / 13 . There 


are three cases, depending on which of the three agents lies 
in between the other twos. 


Case 1. Agent 3ti lies in between xj and X 3 , and without loss 
of generality, we assume that the three coordinates a 02 
and 03 satisfy the condition that 02 < a\ < 03. As both /13 
and /13 are monotonically increasing and have zeros, so if 
we fix agent xa and X3 but move x\ along the a-axis with 
a\ € (02,03), then there is a unique position for agent x\ at 
which it is balanced, i.e, 


fn(dn) ■ (at -a 2 ) = /13W3) ■ («3 -«t) ( 21 ) 

We then define a map £23 : R+ —> R. by requiring 
^ 23 (^ 23 ) • («3 - 02 ) = fn(dn)-{ai-a 2 ) 

= /l 3 (^ 13 ) • (03-a-i) 


and we call £23 the virtual interaction between agents xi and 
if3 induced by /12 and /13. We here note that in this case, the 
map g23 is a function in ■'¥. To see this, we first notice that 
£23 is continuous differentiable, with its derivative given by 


£23 (*3) 

We then notice that 


f[ 2 (dn)f[ 3 (du) 

f[ 2 (dl 2 )+f[ 3 (dl 3 ) 


( 23 ) 


123(^23) =fn(dn) <fu(d 23 ) ( 24 ) 


and hence 

lim / S2i{x)dx < lim / fn(x)dx = — °° ( 25 ) 

d^oJd d->Oj d 

So the map £23 satisfies both conditions Cl and C 2 , and 
hence £23 is contained in &. 

Case 2 . Agent X2 lies in between xi and X3, and we assume 
a 1 < a2 < 03. Similarly, if we fix xx and X3 but move x\ 
along the a-axis with a\ G (—°°,a 2), then there is a unique 
position for agent T] at which it is balanced. Again, we have 


fn(dn) ■ {ai — 02) = f\3(d\3) -(03 — 01) ( 26 ) 


and we define the virtual interaction £23 : R+ —> K by 
requiring 


£23(^23) • (03 — ^2) = fn(dn) ■ (a\ ~02) 
= fu(d 13) ■ (03 ~a\) 

Similarly, the derivative of £23 is given by 


( 27 ) 


£23(^23) 


/T 2 (tfl2)/jf 3 (^13) 
fn(du) + f[3(di3) 


( 28 ) 


but in this case lim,/^ > o<? 23(< : 0 is finite. To see this, we notice 
that both /12 and /13 are strictly increasing, so there is a 
unique do > 0 such that 


Mdo)+fM) = 0 ( 29 ) 


but then, we have 

lim£23(c/) =fn(do) > ( 30 ) 

d->0 

We note that although £23 is not in yet / + £23 will be 
in ■'¥ for any / in &. 


Case 3. Agent x-\ lies in between X\ and X 2 - We then follow 
the same procedure as we did in case 2 to construct the 
virtual interaction £ 23 . Again, in this case £23 is not contained 
in but still we have £13 (t/) > 0 for all d > 0 . 

Equipped with the notion of virtual interaction, we will 
then be able to introduce the reduced system. Let G = (V,E) 
be a TLG of N vertices. Choose a Henneberg construction 
of G, and we assume, along the rest of this paper, that vertex 
1 is the last vertex joining to the graph, via edges ( 1 , 2 ) and 
(1,3) to vertices 2 and 3. Let £23 be the virtual interaction 
between X 2 and X 3 induced by /12 and /13 defined in any of 
the three cases. We then let 

J23(d)-=f23(d)+g23(d) (31) 

and it is clear that /|3 is contained in &. We then define 
a new system of (N — l)-agents by ruling out agent X] , and 
meanwhile replacing /13 with /| 3 . This newly defined system 
will be referred as a reduced system in the rest of this paper. 
As there are three different ways to define £ 23 , and hence 
/| 3 , so there will be three types of reduced systems. 

As later in this paper, we will consider variations of 
interaction laws in both the original system and one of 
its reduced systems. So we find it helpful to introduce a 
map describing the relation between the two ensembles of 
interaction laws. Lor convenience, we let 

f := & 2N ~ 3 

{ a* := & 2N - 5 (32) 

be the collections of ensembles of interaction laws of the 
original system and a reduced system, respectively. We then 
consider maps 

pr-Cl-t ft* (33) 

with i = 1,2,3, each of which is defined by 

P,(/i2,/i3,/23,---) :=C/?3>-) (34) 

The sub-index i of p, indicates in which case the map /,* 3 is 
defined. We here note that each map p, is open, suijective 
and continuous, and we refer readers to the appendix for a 
complete proof of this fact. 

The notion of reduced system will be a useful tool for 
computing the inertia of the Hessian matrix of a critical orbit, 
and for analyzing generic properties associated with the set 
of critical orbits. In particular, it enables us to apply the 
technique of induction as we will see later in this paper. 

VII. The inertia of the Hessian matrix 

ASSOCIATED WITH THE REDUCED SYSTEM 

In this section, we will describe a useful property of 
reduced system which is related to the inertia of the Hessian 
matrix at a critical line configuration. 

We start by introducing some useful definitions and no¬ 
tations. Let p be a fixed line configuration, and we say a 
reduced system is associated with p if the virtual interaction 
£23 is defined with respect to the arrangement of positions 



(42) 


of the three agents x \, x 2 and x 3 in p. Let M be a symmetric 
matrix, and we will let 


n{M) := (n+(M),n 0 (M),n_(M)) (35) 


be the inertia of M. For convenience, we also define a vector¬ 
valued sign function by 

r (1,0,0) if x > 0 

sgn(x) := < (0,1,0) ifx = 0 (36) 

[ (0,0,1) if x < 0 

We will now state the property of the reduced system 
associated with p 

Theorem 2: Suppose p is a critical line configuration, then 
the sub-configuration p*, formed by agents X 2 ,--- ,x„, is a 
critical line configuration of the reduced system associated 
with p. Let <t>* be the associated potential function of the 
reduced system, and let FL p * be the Hessian matrix of <t>* at 
p*, then we have 

n{Hp)-n(H p *) = sgn(-f l2 (d l2 ) - fn(dn)) 
+sgn(-f[ 2 (di 2 )-f[ 3 (d l3 )) 

The expression above we be referred as the inertia formula 
in the rest of this paper. 


It is clear by construction of the virtual interaction that the 
sub-configuration p* is an equilibrium of the reduced system. 
So in the rest of this section, we will focus on the proof of 
inertia formula. We will assume that the line configuration 
p is on the a-axis. By previous computation, we have 


H p = 




(38) 


Let F p * and dF p * be two matrices defined in the same way as 
F p and dF p , yet with respect to the reduced system associated 
with p. Then we have 


H P - = 




(39) 


It then suggests that we relate n{F p ) and n{dF p ) to n(F p *) 
and n(dF p «), respectively. The results are then summarized 
in the statements of lemma [3] and lemma 0] 


Lemma 3: The inertia of F p and F p * are related by n(F p ) — 
n(F p *) = sgn(-f l2 (d l2 ) -fn(du)). 

Proof: Let vj,--- ,v* N _ x be a set of orthonormal eigen¬ 
vectors of Fp* with respect to eigenvalues Ai, • • • . Av 1 . We 
will now use each v* to construct a vector v,- in M. N . Let 
v*- be the j-th entry of vector v?, and as usual we let ap 
be the coordinate of agent xp on the a-axis. We then define 
vt := (a,v*) by adding a scalar a in front of vj\ and the 
scalar a is given by 

a . = (g3-gi)v?i + (fli-g2)va (4Q) 

«3 — a 2 

Notice that in any of the three cases, the virtual interaction 
g 23 satisfies the condition 

£ 23 (^ 23 )' («3 ~a 2 ) = fn(di 2 )-(a\-a 2 ) 

= fn{di 2 ) ■ (<23 -a\) 


This, in particular, implies that 

F p -Vi = Xi-{ 0,0 

for all i = 1, • • • ,N — 1. 

Let e\ be a unit vector with one on its first entry, and zeros 
elsewhere. Let Q be a N-by-N matrix, with e\,vi,--- ,vn-\ 
its column vectors, then Q is a full-rank matrix. We then 
consider a congruence transformation of F p by 

A := Q T F P Q (43) 

By expression ( |42| ), we know that A is a diagonal matrix 
given by 

A = diag(-fi 2 (dn) ~ fi 3 (di 3 ),Xi,- ■ ■ ,Av_i) (44) 


Then by applying the Sylvester’s law of inertia [8], we know 
that the inertia of F„ coincides with the inertia of A. This 
then proves the lemma. ■ 

Lemma 4: The inertia of dF p and dF p * are related by 
n(dF p ) — H(dF p *) = sgn(-f[ 2 {d\ 2 ) -f[ 3 (di 2 )). 


Proof: The proof here will be very similar to the proof 
of lemma |3] Let u \, • • • , iif j be a set of orthonormal eigen¬ 
vectors of dFp* with respect to eigenvalues Ai,• • • ,A^_i. 
Similarly, we let ui := (P,u*) with 


n , = f[ 2 (dn)u* x + f[ 3 (d n )u* 2 
f\ 2 (d\ 2 ) + f[ 3 (d\f) 


(45) 


where u* x and u* 2 are the first and second entries of u* 
respectively. Notice that 


123(^23) 


f[ 2 (dn)f[ 3 (d l3 ) 

f[ 2 (dn) + f[ 3 (dn) 


(46) 


So then, we have 


dFp ■ Ui = A, ■ (0, u*) (47) 

Now let g be a N-by-N matrix, with e\.u \ ,■ ■ ■ , un- i its 
column vectors, then 


Q r dF p Q = diag(-f[ 2 (dn) - f[ 3 (d n ),X h -■ ■, A, V - 1 ) (48) 

The inertia of dF p thus coincides with the inertia of the 
diagonal matrix above. This then proves the lemma. ■ 

Theorem [2] is then established by combining lemma [3] and 
lemma 0 ] 


VIII. Nondegenerate critical orbits of 

LINE CONFIGURATIONS 

Our goal in this section is to prove the next theorem. 

Theorem 5: Let K be the collection of critical orbits of 
line configurations of system 0. then genetically K is a 
finite set and each critical orbit in K is nondegenerate. 

The proof of theorem [5] is divided into two parts. We first 
show that K is generically a finite set, and then we show that 
generically each critical orbit f/ p in K is nondegenerate. 





A. Proof that K is genetically finite 


We will prove this fact by induction on the number of 
agents. 

Base case. Suppose N -- 2, then there is only one critical 
orbit i & p , characterized by the condition that ^ 12 (^ 12 ) = 0 . 
Inductive step. Suppose the lemma holds for N = n with 
n > 2, we then prove for N = n+l. 


We recall that £2* is the collection of ensembles of 
interaction laws in a reduced system. Let £2^ C £2* be defined 
by collecting ensembles of interaction laws by which there 
are only finitely many critical orbits of line configurations. 
Then by induction, the set £2^ contains an open dense subset 
of £ 2 *, and for simplicity, we may assume that £ 2 j itself is 
open and dense in £ 2 *. 


Let £2i be a subset of £2 defined in the same way as £2j in 
£2*, yet with respect to the original system. Let p\, p 2 and 


P 3 be maps from £2 to £2* debited at the end of section VI 
We then consider a set 


3 

£2'i := flPr'^T (49) 

i=\ 


It is clear that £2] contains £2j as a subset. So it suffices to 
show that £ 2 '| is open dense in £ 2 . 

First we show that £2'j is open in £2. Each p^*£2| is open 
in £2 because p, is continuous, and so is their intersection. 
We now show that £2( is dense in £2. Suppose not, there is 
an element ft) in £2 and an open neighborhood U of CO in 
£2 such that I/n£2j = 0 . Then there must be some p, with 
Pi(U) H£2| = 0. On the other hand, the set pi(U) is open in 
£ 2 * because p, is an open map, and £ 2 j is dense in £ 2 *, so 
Pi(U) must intersect £2j which is a contradiction. 


B. Proof that each G p in K is genetically nondegenerate 

It suffices to show that generically hq {H p ) = 3 for all G p 
in K. We again prove this fact by induction on the number 
of agents. 

Base case. Suppose N = 2, then there is only one critical 
orbit G p . Suppose p is on the n-axis, then we have 

F P = (0 0 ) and dF p =f' n {dn)^ (50) 

and hence no(H p ) = no(F p ) +no(dF p ) =2 + 1 = 3. 

Inductive step. Suppose the lemma holds for N = n with 
n > 2, we then prove for N = n+ 1. 

Let £2, be a subset of £2j debned by collecting ensembles 
of interaction laws by which, we have no(H p *) = 3 for each 
critical orbit G p * of line conhgurations. By induction, the 
subset £ 2 , contains an open dense subset of £ 2 *, and we again 
assume that £2^ itself is open and dense in £2*. Similarly, we 
let £22 be a subset of £ 2 i debned in the same way as £ 2 ^ in 
£2*, yet with respect to the original system. Let 

3 

£2t := P| Pr‘^2 (51) 

i= 1 


then by the same reason as £ 2 , in £ 2 , the set £ 2 ', is also open 
dense in £2. Let 

£2? :=£2 2 H£2; (52) 

It then suffices to show that £2" is open and dense in £2(. 

i) Openness. We hrst show that £2" is open in £2. Let ft) 
be an element in £ 2 ", and let @p\ >" ‘ > ^Pk be the critical 
orbits in K associated with ft). Then for each i = l,-- ,k, 
there is an open neighborhood 6 ', of ft) in £ 2 , and an open 
neighborhood V, of G Pi in P such that if ft)' is in Up then there 
is a unique nondegenerate critical orbit Gp in V,. Moreover, 

"i 

the critical orbit vj is an orbit of line configurations. So if 
we let U := nf = 1 f/,, then U is an open subset of £2 contained 
in iff 

ii) Density. We now show that £2, is dense in £2(. Let ft) 
be an element in £2(, we show that any open neighborhood 
U of ft) in £2 intersects £2". Let @p \) ‘" 1 @Pk be the critical 
orbits in K associated with ft). Let p* be the sub-conbguration 
of pi formed by agents j? 2 , ■ ■ ■ ,3c«+i, then each G p * is a 
nondegenerate critical orbit of the reduced system associated 
with p,. By the same arguments as we used to prove the 
openness of £ 2 " in £ 2 , we conclude that there is an open 
neighborhood U of ft) in £2, and an open neighborhoods V, 
of each G Pi in P such that if ft)' is in I/, then there is a unique 
critical orbit of line conhguration G p > in V,. Moreover, each 
Gp* is a nondegenerate critical orbit of the reduced system 
associated with p\. In the rest of the proof, we will bx this 
open neighborhood (/, and we assume that any perturbation 
of ft) is made within U. 

Let W, be a subset of U collecting any ft)' by which 
no(H p i) = 3. It then suffices to show that any open neigh¬ 
borhood (/' of ft) in U intersects W, for each i = 1, • ■ •, k. By 
theorem [ 2 ] we know 

3 < n 0 (H Pi ) < 5 (53) 


and we may assume the worst case where G Pj is degenerate 
with no(H Pi ) = 5, i.e. 


(54) 


f /l2(^12)+/l3(^13) =0 

1 f[ 2 (d l 2)+f[ 3 (d l3 )=0 

We will now perturb ft) to an ft)' in 1+ H U'. For simplicity, 
we only focus on the case where agent x\ is in between x 2 
and X 3 in /?,, but the analysis will be the same in the other 
two cases. Choose a triplet (£12,613,£23) in IIjLjC 1 (K+,K) 
such that they satisfy the condition 

( £12(^12) = £13(^13) = —£23(^23) > 0 
\ e[ 2 {d n ) + e’ n iydi 3 )>Q 

and we assume that (£12,£13,£23) are small enough such that 
ft)':= ft)+ (£i 2 ,£ 13 ,£ 23 ,0, •■•,()) (56) 


(55) 


is contained in I/'. Then by construction, the perturbed 
critical orbit Gp coincides with the original one G Pi . But 
by this perturbation, we have 

(fn(d l2 ) + £12(^12)) + {fri(dif) + £13(^13)) 

= £i2(ft f i2)A(i2 +£13(^13)7^13 > 0 



and also 


( f ' nid ^ + e ' ntdn )) + (/13 (^13)+ £13(^13)) 
= e( 2 (^i 2 ) + £13(^13) > 0 


suffices to show that for any function /| 3 in &, there is a 
triplet (/i2,/i3,/23) in such that 

f23=f23+823 (62) 


So then no(H p /) = 3 after perturbation, and hence (o' is 
contained in Wj. This then completes the proof. 

IX. Proof of the main theorem 

Let G = (V,E) be a TLG of N vertices, and let S& be 
the collection of all triangulated Laman subgraphs of G. 
Let G' = iV',E') be an element in Sg, and we consider a 
sub-system induced by G', i.e, the sub-system formed by 
agents x,- with i £ V' , and {fjk\(j,k) £ E'} is the ensemble 
of interaction laws. Let 

■= & |£1 (59) 

be the collection of ensembles of interaction laws associated 
with the sub-system. Let fig/ C fig/ be the subset collecting 
any ensemble of interaction laws by which the associated 
potential function of the sub-system is an equivariant Morse 
function. By theorem [ 5 ] the set fig/ is open dense in fig/. 
Now for convenience, we let 


1 

fej. 

II 

O 

l 

a 

(60) 

and we define a subset of fl by 


fl := P| (fig/ x fl_g/) 

(61) 




with g23 the virtual interaction induced by /12 and /13. Let 

Md):=\fhV) 

fn{d)=f n {d):=\f^{2d) 

for all d > 0 . We then verify that all /12, /13 and /13 are in 
and £23 = fi 3/2 is the virtual interaction induced by /12 
and /i 3 . Thus, the expression ([ 62 } is satisfied by our choice 
of /12. /i3 and / 23 . 

2 . Proof that p\ is open 

Let U be an open set in fl, let U* ■= p\(U) and we show 
that U* is an open set in fl*. Let < 0 * be in TJ*, and let 
&> = (••• ,fij ,-•■) be in U with pi(ft)) = ft)*. Choose open 
neighborhoods Vjj of fj in & such that Yl( i ^ eE V i j is in U. 
Let 

E* :=,E — {( 1 , 2 ),( 1 , 3 )} ( 64 ) 

and let U\ be a subset of U defined by 

U\ := /12 x /13 x LI(j J ) e £»V'23 ( 65 ) 

Let U* := pi(f/i), then we have 

U\ :=H^j) eE *Vij - H^p eE *fij + < 0 * (66) 

So t/j* is an open neighborhood of ft)* contained in U*. 


Since fl is a finite intersection of open dense subsets of fl, 
and hence fl itself is open dense in fl. 

It now suffices to show that for each element ft) in fl, 
the associated potential function is an equivariant Morse 
function. Let be a critical orbit, and let {pi}f = \ and 
{G,}^i be sub-configurations of p and sub-graphs of G 
respectively, associated with the canonical partition. Then 
each 0pt is a critical orbit of line configurations of the 
sub-system induced by G,-. As each G; is a TLG, so by 
construction of fl, we have hq{H P i ) = 3 for all / = 1, • • • ,m. 
We then apply the Morse-Bott index formula to conclude that 
no(H p ) = 3. In other words, we have just showed that each 
critical orbit fj p is nondegenerate. As there are only finitely 
many critical orbits, so the associated potential function is 
an equivariant Morse function. 

Appendix 


3. Proof that p 1 is continuous 

Let 77 : JF 2 —>• & be the map defined by sending a pair 
C/ 12 ,/ 13 ) to the virtual interaction £ 23 , it then suffices to 
show that 77 is continuous. The proof will be divided into 
two parts, we will first fix a distance d > 0 , and prove that 
the map T),i : cF 1 -O R 2 defined by 

rid :(fn,fi 3 )^ (823(d),g' 23 (d)) (67) 

is continuous. The continuity of 77 will then be established 
afterwards. 

i) Continuity of r\ t i 

Let (/i 2 ,/i 3 ) be a pair contained in ^ 2 , and let (dn,d\i) 
be the unique pair of positive numbers with 

d\2 + dn = d 

( 68 ) 

fn{dn) = fn(d\f) 


We here will prove that each map p ( , with i = 1,2,3, 
defined at the end of section VI is open, surjective and 
continuous. For simplicity, we will only focus on the proof 
about p 1 , the same arguments can be applied to the other 
two cases. 


1. Proof that p\ is surjective 

Let ft)* be an element in fl*, and we now show that there 
is an element ft) in fl with p(ft)) = ft)*. By construction, it 


It is clear that for fixed d, the pair (^ 12 ,^ 13 ) is a function of 
C/12,/13) only. Notice that 


823(d) = fn(dn) 

1 f' l2 {dn)+f' u {<hi) 


(69) 


So it suffices to show that the pair (c/12,^13) is continuous 
in (/12, / 13). We shall use the inverse function theorem to 
establish this fact. 




As continuity is a local problem, so it suffices to consider 
variations of /12 and /13 over some neighborhoods of dn 
and d\ 2 . Let 7 be any closed neighborhood of 0 in R, and 
let /12 and /13 be closed intervals defined by 


hi :=dn + I 
hi -=d 13 —7 


(70) 


So then 7i2 and 7x3 are closed neighborhoods of d\2 and d\ 3 
respectively. We now consider a C-map 

E, :7xC 1 (7 12 ,M) xC 1 ^,®) ->• R (71) 


defined by 

<i;(x,/zi2,/ii3) :=hn{dn+x)-hn(di-i-x) (72) 


Let /121/12 and /i3|/ 13 be restrictions of /12 and /13 to 7x2 and 
7 j 3 respectively, then we have 


^(/l 2 |/ 12 ,/l 3 |/ 13 , 0 ) = 0 

= /( 2 (^ 12 )+/( 3 (^ 13 ) > 0 


K 

dx 


(73) 


(/l2l/ 12 >/l3l/ 13 i0) 


So by the inverse function theorem for the Banach spaces, 
we know that there exist 


1) . an open neighborhood V12 of fn\i l2 in C 1 (7i2,R); 

2) . an open neighborhood Vi3 of /131/ 13 in C^^R); 

3 ) . an open interval (—t,t), with t > 0 , in R; 

4 ) . a unique C'-map a : Vn x V13 (~t,t) 


such that 

£(/U2,/L3,C7(/zi2,/zi3))=0 (74) 


for all (/z 12, 7 z13) in V12 x V13. This then shows the continuity 
Of (c/12,^ 13 ) at (/l 2 ,/l 3 )- 
ii) Continuity oft J 

Let (/i2,/i3) be any pair in d? 2 , we will prove the 
continuity of 77 at (/i2,/t3)- As the distance d is no longer 
fixed, so to avoid confusion, we will write the argument d 
explicitly out, e.g, we will write d\ 2 (d) and d\ 2 (d) instead 
of d 12 and c/13. 

We start by defining a family of closed intervals 1 (d), 
parametrized by d over R + . Let 


I(d):=[~a(d),p(d)} (75) 

with both a and functions in C°(R + ), and we let 

( hi(d) :=d l2 (d)+I(d) 

\ hi(d):=d n (d)-I(d) 

be closed neighborhoods of d\ 2 (d) and d 13(d) respectively. 
Now for each d > 0 , we define two subsets of R + by 


J\i(d ){ d' £ R+|c/ £ l n (d')} 
J 13(c/) := {d f £ R + |c/ £ 7i3(c/ / )} 


( 77 ) 


for all d > 0 . 

We will now show that for any 8 in C°(R+), there exists 
an e in C°(R + ) such that 

Be(fn) xB e (f 13 ) C 1 T 1 CBsfe)) (78) 

First by the proof of continuity of rj d , we know that for a 
fixed d > 0 , there exists a number rj > 0 such that if 

max{|e(x)| + |e'(x)||xe7i2(c/)U7i 3 (c/)} <r d (79) 

then for any (hn,h\l) in B e (f n ) xB £ (f l3 ), we have 

|t] ( i(/zi2,/zi3) —I?d(/l2,/l3)|i < 8(d) (80) 

where | • |i denotes the 1-norm of a vector in M 2 . Let 

r d :=min{r d i\d' £ J\ 2 (d)CJ 13(d)} (81) 

As J\ 2 (d) and J\i(d) are compact, so then r d > 0. Now we 
choose an e in C°(R+) such that e(d) < r d for all d > 0. 
This then implies that 

t](h\ 2 ,hi3) £B 5 (g 2 3) (82) 

for all (Zzi 2 ,/zi 3 ) in B e (f n ) x B e (/i 3 ). 
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We then assume that a and /j are chosen such that 

1. both In(d) and I\3(d) are in R + ; 

2 . both J\ 2 (d) and J\i(d) are compact 



